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Solar neutrino fluxes depend both on astrophysical and on nuclear physics inputs, namely on the 
cross sections of the reactions responsible for neutrino production inside the Solar core. While the 
flux of solar *B neutrinos has been recently measured at Superkamiokande with a 3.5% uncertainty 
and a precise measurement of ^Be neutrino flux is foreseen in the next future, the predicted fluxes 
are still affected by larger errors. The largest nuclear physics uncertainty to determine the fluxes 
of *B and ^Be neutrinos comes from the '^He(a,7)^Be reaction. The uncertainty on its S-factor is 
due to an average discrepancy in results obtained using two different experimental approaches: the 
detection of the delayed 7 rays from ^Be decay and the measurement of the prompt 7 emission. Here 
we report on a new high precision experiment performed with both techniques at the same time. 
Thanks to the low background conditions of the Gran Sasso LUNA accelerator facility, the cross 
section has been measured at Ecm = 170, 106 and 93 keV, the latter being the lowest interaction 
energy ever reached. The S-factors from the two methods do not show any discrepancy within the 
experimental errors. An extrapolated S(0)= 0.560±0.017 keV barn is obtained. Moreover, branching 
ratios between the two prompt 7-transitions have been measured with 5-8% accuracy. 

PACS numbers: 25.55.-e, 26.20.+f, 26.65.+t 

Keywords: ^He(a, 7)^^66, p-p chain, direct measurement, underground accelerator 



Forty years ago, John Bahcall and Raymond Davis 
started to explore the solar interior by studying the neu- 
trinos emitted by the Sun The results of the first 
neutrino detection experiment Q originated the so called 
solar neutrino puzzle, consisting in a deficit of measured 
neutrinos with respect to the theoretical predictions of 
the Standard Solar Model (SSMV After thirty years of 
experiments, SNO and Kamland 0, Q observed neutrino 
oscillations and proved that the missing solar electron 
neutrinos actually change their flavour during the travel 
to the Earth. This closed the neutrino puzzle. There- 
fore, the high precision measurement of ^B neutrino flux 
, together with the foreseen measurement of ^Bc neu- 
trinos 0, can now be used to understand physical and 
chemical properties of the Sun, provided that nuclear rc- 
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action cross sections are known with similar accuracy [7| . 
The ■^He(Q!,7) ""Be reaction is the onset of the ^Be and ®B 
branches of the proton-proton (p-p) chain of hydrogen 
burning. The 9% error [8| on its cross section is presently 
the main nuclear physics uncertainty on the prediction of 
^Be and ^B neutrino fluxes 0]. 

At stellar energies the 3He(Q;,7)'^Bc cross section 
(t(E) drops exponentially with the energy and can be 
parametrized as: 

E 

where S(E) is the astrophysical factor, 77 is the Sommer- 
feld parameter p^ . and E is the center of mass energy. 

3He(a,7)^Bc is a radiative capture reaction (Q-value: 
1.586 MeV) into the first excited state (E^=429 keV) 
and the ground state of ^Be that subsequently decays 
by electron capture into ^Li with a terrestrial half life of 
53.22±0.06 days 

In the last forty years the ^He(a,7)''Be reaction has 
been measured using two techniques. In the first ap- 
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proach direct a-capture 7-rays were detected (prompt 7 
method) [H, [H, [H, O [H , while, in the second, 
the delayed ''Be-decay 7 rays were counted (activation 
method) [T^. [l9l. [20l. w\ . Previous activation results are, 
on the average, 13% higher than prompt 7 data and this 
is the origin of the large uncertainty quoted on the re- 
action cross section Q. Up to now, no explanation has 
been obtained for this discrepancy that could be due ci- 
ther to systematic experimental errors (angular distribu- 
tion, branching ratio eflFects, parasitic reactions produc- 
ing ""Be) or to the existence of a non radiative capture 
(EO monopole) [IJ . Recently, the discrepancy on the ex- 
trapolated S(0) has been reduced by an activation study 
at 420<E<950 keV ^ to 9%. High accuracy (4%) ac- 
tivation data were obtained also at LUNA, at center of 
mass energies down to 106 keV [2^ [l^j- Nevertheless, 
high accuracy prompt gamma data are also needed to 
verify the claimed discrepancy. 

We have performed a new high accuracy measure- 
ment using simultaneously prompt and activation meth- 
ods with the same experimental setup. The experiment 
has been carried out using the underground LUNA 400 
kV accelerator [2^ at the Gran Sasso National Labora- 
tory (LNGS). 

Three couples of cross section values have been mea- 
sured (prompt 7 and activation) at Eq=220, 250 and 400 
keV. A sketch of the interaction chamber is given in Fig. 
[TJ The a beam enters the '^He extended windowless gas 
target [23 | through a 7 mm diameter collimator and is 
stopped on a detachable copper disk that serves as the 
primary catcher for the produced ''Be and as the hot side 
of a calorimeter [2^. The latter measures the beam in- 
tensity (about 250 fj,A) with an accuracy of 1.5%. The 
high beam current decreases the "^He density along the 
beam path this effect has been monitored with a 

silicon detector by double Rutherford scattering provid- 
ing an accuracy of 1.3% on the gas density determination 
[28l |. The same detector is also used to measure the gas 
contamination (mainly N2) that has remained below (2.7 
± 0.3)%. 

Prompt 7 rays are counted with a 135% ultra low- 
background HPGe detector shielded with 5 cm of OFHC 
copper and 25 cm of lead. The detector and the shield are 
enclosed in a sealed plastic box flushed with dry N2 to 
reduce ^^^Rn background. Thanks to the underground 
environment where cosmic muons are strongly reduced 
[2^, the shielding suppression factor is of five orders of 
magnitude for 7 rays below 2 MeV. A lead collimator is 
positioned inside the target chamber to collect mostly 7 
rays emitted at 55°. At this angle the contribution of 
the second Legcndrc polynomial in the angular distribu- 
tion expression, vanishes. Therefore the inner collimator 
reduces the systematic error due to prompt 7 angular 
distribution uncertainties and also shields the detector 
from possible beam-induced radiation coming from the 
entrance collimator and the calorimeter cap. The effec- 
tive target length seen by the HPGe detector is approx- 
imately 12 cm, corresponding to an energy loss AE = 3 
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FIG. 1; Schematic view of the interaction chamber with the 
position of the HPGe detector and of the 100 ^m silicon detec- 
tor used for ^He density monitoring. The distance between 
the entrance collimator and the calorimeter is 35 cm. The 
thickness of the internal collimator is 3 cm for the Lead part 
and 1.6 cm for the Tungsten part. 



keV at Ptarget=0.7 mbar (value used in all the runs) and 
E„=400 keV. 

The photopeak detection efficiency is determined by a 
Monte Carlo code ^ calibrated with ^"Co and ^^^Cs ra- 
dioactive point-like sources moved along the beam path. 
The Monte Carlo reproduces the experimental cfhcicncy 
within the source activity uncertainties (1.5%). The 
spectra collected at Ea= 220, 250 and 400 keV, with a to- 
tal charge of 637, 407 and 113 C respectively, are shown 
in Fig. [2] together with laboratory background. Beam 
induced gamma-ray background has been measured with 
^He gas in the target at Eq= 400 keV: no difference with 
laboratory background has been observed. In the data 
analysis only the two primary transitions at E.y=Q-|-Ecm 
and E-,=Q+Ec,n— 429 keV, have been considered. The- 
oretical angular distribution functions are calculated by 
[3H down to 210 kcV interaction energy. A linear ex- 
trapolation of the curves of has been done and the 
coefficients of the Lcgendre polynomials adopted in the 
detection efficiency calculation are: ai= - 0.05 and ai = 
for the transition to the ground and to the first ex- 
cited state, respectively and a.2= - 0.1 for both transi- 
tions. These values are in agreement with recent theo- 
retical predictions 3^. With the Monte Carlo code, we 
have conservatively varied 100% both ai and 3.2 coeffi- 
cients obtaining a global 2.5% variation of the detection 
efficiency. This value has been assumed as a systematic 
uncertainty and turns out to be the major contribution 
to the error budget of the prompt 7 method. The in- 
beam runs provide accurate branching ratios between the 
two transitions cr(DC^ 429)/(t(DC->0): 0.417 ±0.020, 
0.415±0.029 and 0.38±0.03 at E„ = 400, 250 and 220 
keV, respectively. Our values are consistent with, but 
moreprecise than, latest branching ratio measurements 
[Tq|and are in agreeement with theoretical calcula- 
tions I33II. 
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FIG. 2: 7-ray spectrum at E„ = 220, 250 and 400 keV com- 
pared with natural laboratory background (in grey) normal- 
ized to beam-measurement live times (respectively: 31.2 days, 
21.3 days and 4.8 days). Arrows indicate the primary transi- 
tion peaks to the first excited state and to the ground state. 

During the runs in which prompt 7-rays are detected, 
the ''Be nuclei produced inside the gas target get im- 
planted into the removable calorimeter cap. After each 
run, the cap is dismounted and moved to LNGS under- 
ground low-activity counting facility [2^ . Details and ac- 
curacy of the activation method at LUNA arc discussed 
elsewhere [1^ HI]. Since we have simultaneously used 
the same beam and target for both methods, some sys- 
tematic uncertainties (beam intensity, target density and 
purity) cancel out in the comparison between the two 
techniques. 

Results are reported in Table U and shown in Fig. [3] 
together with all previous literature data. For each cou- 
ple of data (prompt 7 and activation) obtained at the 
same Ea, the effective energy (Eg//), calculated as de- 
scribed in [2^, is slightly different. Indeed the target of 
the prompt 7 experiment, defined by the inner collima- 
tor (Fig. [1]), is a fraction of the whole target contributing 
to ^Be production. The Ee// difference corresponds to 
an S-factor change smaller than 0.1% according to the 



energy dependence given in [3J|. In the comparison be- 
tween prompt and activation S-factors, we have therefore 
neglected the Ee/ / differences and considered a total un- 
certainty given by the statistical and reduced systematic 
errors summed in quadrature (Table [l|. The mean per- 
centage difference between the S-factor values in Table |T] 
(AS=(Sa-Sp)/((Sa+Sp)/2)) is AS„= -0.014±0.042. This 
result limits to +2.8% (maximum AS™ value at Itr level) 
possible non-radiative contributions to the reaction cross 
section. S-factor activation values at Eq=400 and 250 
keV arc compatible with those previously obtained at 
LUNA with the same setup [2^ [2J| . Considering the av- 
erage of the new and old activation values at the same 
beam energy, AS„i does not change. A simultaneous 
measurement with both activation and prompt 7 tech- 
nique at energies around Ecm=l McV, where the oldest 
activation experiments p^.[l9j were performed, would be 
useful to look for non-radiative contributions in a higher 
energy region than the one explored in the present exper- 
iment. To deduce the extrapolated S(0), the fit of has 
been rescaled using the present activation and prompt 7 
data separately. The weighted average between the two 
S(0) values has been calculated adopting as weights the 
statistical error obtained from the fit and the reduced 
systematic error. We get S(0)=0.560±0.017 kcV barn 
where the final uncertainty also includes the systematic 
error common to the two methods. Performing the same 
calculation considering also the most recent and very ac- 
curate results from [21I, [23| or using the theoretical 
function [s^ adopted in the NACRE compilation [sB], 
instead of the R-matrix fit by the extrapolated S(0) 
changes less than 1%. Low energy accurate data in fact 
minimize the uncertainty upon extrapolation. However, 
a refined measurement of the slope of the S-factor in a 
wide energy range would be useful to confirm theoretical 
calculations reducing the uncertainty on the extrapolated 
S(0). The uncertainty on the predicted ®B neutrino flux 
due to S34 is now reduced from 7.5% to 2.4% and the to- 
tal uncertainty, including astrophysical parameters, goes 
from 12% to 10% [s^]. Similarly, the uncertainty on ^Bc 
predicted flux goes from 9.4% to 5.5%, being the contri- 
bution of S34 error reduced from 8% to 2.5% (37j . 
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FIG. 3: Overview of all available S-factor values for the ■^He(a,7)^Be reaction. Filled and open circles: present work. Prompt-7 
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[1^ . Activation data: filled diamonds [l9l |. crosses . filled squares [2l|, filled triangles [l^, [Ij. Dashed line: most recent 
R-matrix fit [s^l, solid line: fit [sj] normalized to present data. In the inset a zoom of prompt 7 (filled circles) and activation 
(open circles) data obtained in this work. 
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Ea 


Method 


Ee/7 




S(Ee//) 


AS Stat. 


AS syst. 


AS red. syst. 


(keV) 




(keV) 


(nbarn) 


(keV barn) 


(keV barn) 


(keV barn) 


(keV barn) 


400 


P 


170.1 


10.25 


0.510 


0.008 


0.019 


0.015 


400 




169.5 


10.00 


0.507 


0.010 


0.015 


0.010 


250 


P 


106.1 


0.588 


0.518 


0.014 


0.019 


0.016 


250 


a 


105.7 


0.546 


0.493 


0.015 


0.015 


0.011 


220 


P 


93.3 


0.235 


0.527 


0.018 


0.020 


0.016 


220 


a 


92.9 


0.232 


0.534 


0.016 


0.017 


0.013 



TABLE I: Cross section and S-factor results with corresponding uncertainties for the prompt (p) and activation (a) methods. 
In the last column the reduced systematic uncertainty is reported where contributions common to the two methods have been 
subtracted. 



